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ABSTRACT 

The influence of primary orientation on the elastic 
response of a [001]-oriented nickel-base single-crystal 
superalloy, PWA 1480, was investigated under 
mechanical, thermal, and combined thermal and 
mechanical loading conditions using finite element 
techniques. Elastic stress analyses were performed using 
the MARC finite element code on a square plate of 
PWA 1480 material. Primary orientation of the single 
crystal superalloy was varied in increments of 2°, from 
0° to 10°, from the [001] direction. Two secondary 
orientations (0° and 45° ) were considered, with respect 
to the global coordinate system, as the primary 
orientation angle was varied. The stresses developed 
within the single crystal plate were determined for each 
loading condition. In this paper, the influence of the 
angular offset between the primary crystal orientation 
and the loading direction on the elastic stress response 
of the PWA 1480 plate is presented for different loading 
conditions. The influence of primary orientation angle, 
when constrained between the bounds considered, was 
not found to be as significant as the influence of the 
secondary orientation angle, which is not typically 
controlled. 


INTRODUCTION 

Nickel base superalloys have attracted considerable 
interest for use in rocket and gas turbine engines due to 


the superiority of their material properties at high 
temperature operating conditions (Chandler, 1985). 
Initial development of these alloys began in the 1960’s 
for application as turbine blade and vane materials. 
Pratt and Whitney Aircraft was among the first to 
manufacture single crystal blades, using PWA 1409 
which is the single crystal version of Mar-M200 (Kear 
and Piearcey, 1967). Since that time many other engine 
companies have developed similar alloys. However, the 
alloy of interest in this work is the single crystal 
superalloy PWA 1480 (SC PWA 1480) developed by 
Pratt and Whitney. This superalloy, which is a well 
characterized nickel-base cubic single crystal, was 
selected by Pratt and Whitney as the blade material for 
the Alternate Turbopump Development Program for the 
Space Shuttle Main Engine (SSME). 

The turbine blades of SC PWA 1480 are directionally 
solidified such that the span direction is along the low 
modulus [001] crystallographic direction to enhance 
creep and thermal fatigue resistance. There is some 
variation in the [001] orientation of the single crystal 
from one blade to another due to the manufacturing 
process. However, it is generally controlled to be within 
10° of span direction for all of the turbine blades. The 
transverse orientation or the secondary crystal 
orientation [010] is not controlled, and is randomly 
oriented with respect to a fixed geometric axis in the 
turbine blade. However, using a seed crystal during 
solidification will allow control of the secondary 
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orientation (Duhl, 1989). The single crystal superalloys 
exhibit cubic symmetry. Three independent elastic 
constants are required to describe the elastic behavior of 
a cubic single crystal (Nye, 1957 and Yang, 1984). Their 
elastic stress-strain relationships are highly orientation 
dependent (Bowen, et al., 1986 and Abdul-Aziz, et al., 
1989). In a previous study, the influence of the 
secondary orientation on the elastic response of SC 
PWA 1480 was documented (Kalluri, et al., 1991). 

The objective of this work is to determine the 
influence of the primary orientation on the elastic 
response of [OOlj-oriented, SC PWA 1480 under 
thermal, mechanical, and combined thermal and 
mechanical loads by conducting finite element stress 
analyses. The analyses were conducted by varying the 
primary orientation [001] in increments of 2° from 0® 
to 10° with respect to a global coordinate system. Two 
secondary orientations of 0° and 45° with respect to the 
global coordinate system were found to represent two 
extreme conditions in a previous study (Kalluri, et al., 
1991), and were selected for investigation in the present 
study. The stresses developed within the SC PWA 1480 
material were determined for each loading condition 
and the influence of the angular offset between the 
primary crystal orientation and the global coordinate 
system is discussed for the two secondary orientations 
considered. 


PROBLEM DESCRIPTION 

The turbine blade of a rocket engine turbopump is 
schematically represented in Fig. 1. This hollow blade is 
made of SC PWA 1480 and is directionally cast such 
that the [001] crystallographic direction (primary 
orientation) is offset with respect to the Z-axis by an 
angle no greater than 10°. In this design, the hollow 
core extends below the platform into the shank region. 
In Figure 1, XYZ is the global coordinate system, with 
the Z-axis along the span, Y-axis along the chord and 
the X-axis along the thickness of the blade. 

The turbopump blades in a cryogenic liquid propellant 
rocket engine are cooled except for their airfoils which 
are typically exposed to the hot gas flow path that can 
induce high thermal gradients both along the spa n an d 
through the thickness. For this hollow core blade, the 
thermal gradient through the thickness may be of the 
order of 875° C/cm, while being much milder along the 



FIG. 1: SCHEMATICS OF SSME TURBINE BLADE 
AND THE SQUARE PLATE USED IN THE 
ANALYSIS. 


span. These gradients are quite high compared to 
aircraft turbine engine blades. For example, the first 
stage high pressure fuel turbopump blade of the SSME 
is exposed to an extremely high temperature 
environment where the gas temperature reaches 
approximately 2000 “C at start-up conditions. Lee et al. 
(1991) have conducted a 3-D elastic finite element 
analysis to investigate the stress state of a single crystal 
hollow core blade with respect to the secondary 
orientation variations. The temperature bounds used in 
the current analysis are in agreement with the 
temperature range used in the investigation by Lee, et 
al. (1991). 

The turbopump blade has complex geometry which 
makes the thermal and mechanical stress analyses under 
these severe operating conditions very complicated. 
Therefore, to investigate the effects of both the primary 
and secondary orientations under thermal and 
mechanical operating conditions, an analysis of a square 
plate made of SC PWA 1480, which is subjected to 
similar loading conditions will be employed (Fig. 1). A 
plate with the dimensions of 25.4 mm by 25.4 mm by 3.2 
mm was selected for the parametric finite element stress 
analyses. 
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(a) ROTATION ABOUT THE Z-AXIS 


systems represent the orientations of the primary and 
secondary crystallographic directions, respectively. The 
matrices of direction cosines that relate the global 
coordinate system (XYZ) to the crystal coordinate 
system (X c Y c Z c ) (Fig. 2) are given in the following 
relationships: 

a) Rotating about the Z c -axis 
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The transformation matrix is simply obtained by 
multiplying the above two matrices, Equations (1) and 
(2). The result of this multiplication is represented in 
the following: 


(b) ROTATION ABOUT THE X'-AXIS 

FIG. 2: CRYSTAL AND GLOBAL COORDINATE 
SYSTEMS 
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The choice of the plate dimensions was based on 
convenience and simplicity, except for the thickness 
which represented the wall thickness of the hollow core 
blade’s airfoil-shank region. The crystal (X^pZ^ and 
global (XYZ) coordinate systems for the plate are 
shown in Fig. 2 The crystal coordinate system X c Y c 
Z c is first rotated about the Z c -axis by an angle 8 to 
obtain the X' Y' Z ' coordinate system. The X' Y' Z' 
coordinate system is then rotated about the X'-axis by 
an angle * to obtain the global coordinate system, XYZ 
The rotation angles * and 6 among the coordinate 


With, 

P n = COS0 

P 12 = sin0 
Pl3 = 0 
P 21 = -cos*sin 0 

P 22 = COS*COS6 

P23 = sin* 

P 31 = sin*sin 0 
P 32 = -sin*cos 0 
P33 = cos* 


( 3 ) 
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Equation (3) is the direction cosine matrix and (Jy 
(i,j=13) are the corresponding direction cosines. 

ANALYTICAL PROCEDURE 


A. Elastic Stiffness Coefficients 

To describe the elastic behavior of a cubic single 
crystal such as the SC PWA 1480 within the 
crystallographic coordinate system, three independent 
elastic stiffness coefficients (C n , C 12 and C^) are 
required (Nye, 1957). This stiffness matrix for a cubic 
single crystal in the crystal coordinate system is given in 
Equation (4). 
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where [O'] is a 6x6 stiffness matrix in the global co- 
ordinate system. The stiffness coefficients for a given set 
of the primary and secondary orientation angles form 
a 6x6 matrix that is defined in terms of the directional 
cosines (Eq. (3)) and is documented by Leiberman and 
Zirinsky (1956). Thus, Equations (4) and (5) define the 
stiffness matrix for a cubic single crystal in a global 
coordinate system (XYZ) that is offset by two angles, 
the secondary orientation angle 0 and the primary 
orientation angle ^r, to the crystal coordinate system 
(X.V C Z C ; Fig. 2). The relationships among the stresses 
and strains in the global coordinate system with the 
transformed stiffness matrix are given by: 

M - P'tt (6) 


where, 

= ^ a xx °yy °zz *yz X zx x xy ^ 


and, 

t e = { y e zz fyz Tec Y»y ^ 


The elastic behavior in the global coordinate system 
can be described by transformation of the stiffness 
tensor from the crystal coordinate system to the global 
coordinate system (Fig. 2) according to the law of 
transformation of fourth rank tensors. Liebermen and 
Zirinsky (1956) developed a method for transforming 
the stiffness coefficients of a single crystal between two 
coordinate systems that are related by a matrix of 
direction cosines, Equation (3). Their method was 
applied to derive the elastic stiffness coefficients of a 
cubic single crystal in the global coordinate system. The 
stiffness matrix of the cubic single crystal within the 
global coordinate system [C'] is given by the following 
relationship, (Eq. (5)): 


[C i = m[Qm T 


(5) 


In Equation (7), o ff a u and t are the 

engineering normal and shear stresses, and in Equation 
(8), tyy, c a and y^, y^ y^ are the engineering 
normal and shear strains, respectively. The three 
independent elastic stiffness coefficients for the SC 
PWA 1480 material at various temperatures were 
obtained from Kalluri et al. (1991). 

The secondary orientation angles, of 0® and 45®, were 
selected for investigation because a previous study with 
the same exact boundary conditions conducted on SC 
PWA 1480 (Kalluri et al, 1991) showed that the 
stiffness coefficients have a periodicity of 90®, with the 
secondary orientation angle, 0. It was also found in that 
studv that the stresses generated were maximum at 


0=45® and minimum at 0=0®, which made these two 
secondary orientations to be the logical choices for 
future investigations. However, for the primary 
orientation angle, i|r, the range chosen was from 0® to 
10® because manufacturing tolerance is typically 10® 


for ijr. 
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B. FINITE ELEMENT ANALYSIS 

The elastic stress analyses were carried out using the 
MARC Finite Element code. These calculations were 
initiated by setting up a physical representation of the 
problem through the construction of a finite element 
model for the square plate geometry. The model 
consisted of 500 isoparametric eight-noded solid 
elements with 1040 nodes and was generated using 
PATRAN (1989). Ten elements were generated along 
the sides and five elements were used through the 
thickness of the plate, (Fig. 3). 



FIG. 3: FINITE ELEMENT MESH OF THE 
SQUARE PLATE. 


The stress analyses were conducted by running a heat 
transfer analysis followed by a subsequent elastic stress 
analysis. In order to ensure equilibrium and avoid rigid 
body motion, the following mechanical boundary 
conditions were imposed (Fig. 3): (1) the node at D was 
fixed in the X, Y, & Z directions, (2) nodes along the 
edge AD were fixed in the Y & Z directions, (3) nodes 
along the edge DC were fixed in the X & Z directions 
and (4) nodes along edges AB and BC were fixed in the 
Z direction. Enforcement of all these boundary 
conditions was implemented by invoking the Fixed 
Boundary option in the finite element program. Five 


different loading cases were considered in this work: 
mechanical, thermal (two cases), and combined 
mechanical and thermal (two cases). The previously 
listed boundary conditions were employed in all five 
loading cases. 


1. Heat Transfer Analysis. Thermal loads with 
gradients along both X and Z directions were applied to 
the model. Figures 4(b) and 4(c). Choice of imposing 
the thermal gradients along these two directions was an 
attempt to simulate gradients through the thickness and 
along the span of a turbine blade. The other loading 
case involved a distributed mechanical load of 137.9 
MPa along the Z-direction, Figure 4(a). Other loading 
conditions investigated were combined mechanical and 
thermal loads for the two directions considered (X and 
Z). 

The single crystal PWA 1480 is a thermally isotropic 
materiaL All physical properties used for the analyses 
including the mean coefficient of thermal expansion, 
thermal conductivity, and specific heat were documented 
by Kalluri et aL (1991). The thermal boundary 
conditions imposed were representative of those 
typically experienced by rocket engine turbine blades. 
Figure 4 shows the thermal gradients employed, which 
closely represent the SSME turbine blade temperatures 
at full power level conditions. Plate surfaces, noted on 
Figure 4(b) and (c) with the designated temperatures 
assigned to them, were added to the thermal model in 
MARC Finite Element code and the rest of the faces 
were assumed to be insulated. Output from the thermal 
analysis was used as an input for the subsequent elastic 
stress analyses. 


2. Stress Analysis. The stiffness coefficient matrices 
for the various primary orientations (<r=0° to 10° in 
increments of 2 s ) and for the selected secondary 
orientations (0=0® and 45°) were all incorporated in 
the finite element calculations through the MARC 
Finite Element code. In addition, the temperature 
dependencies of all the physical properties were used. 

To maintain consistency and to compare with the study 
conducted previously on secondary orientation effects 
(Kalluri et aL, 1991), the elements shown in Fig. 3 were 
selected to document the variations in the stress 
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(c) THERMAL GRADIENT ALONG X-AXIS 


FIG. 4: LOADING CASES IMPOSED ON 
THE SQUARE PLATE. 


components as the primary and secondary orientations 
were varied. The reason behind the choice of these 
elements, as described by Kalluri et aL (1991), was 
mainly to capture the locations that represent the true 
effects of the loading imposed, and not the effects of 
thermal and mechanical constraints applied. These 
elements are: Element 295, located near surface EFGH 
for the thermal load in the Z-direction case, where a 
thermal gradient is imposed through faces ABCD and 
EFGH, element 445 for the thermal load in the X- 
direction case (thermal gradient through faces DCGH 
and ABFE) and element 295 for the mechanical loading 
case. 


DISCUSSION OF ANALYTICAL RESULTS 

In this study, the influence of both primary and 
secondary orientation angles on the elastic response of 
SC PWA 1480 was determined under different thermal 
and mechanical loading conditions. For each element 
and loading case, the individual stress components at 
given primary and secondary orientation angles f and 6 
were normalized with the corresponding stress 
components obtained for 4r=0 o , 0=0® or 45®, i.e. 

( a xx)*y(°xx)o,e> ( 0 yy)*y( 0 yy)o,e> ( T xy)* V( T xy)o,6 > and 50 
on. For each secondary orientation angle considered, 

the individual as well as the normalized components 
were plotted against the primary orientation angle for 

each element and loading condition. Figures 5 to 10 
show these plots as a function of the primary 
orientation angle i(r, for fixed secondary orientation 
angles of 0=0® and 45®. 

Figure 5 1 shows the influence of the primary 
orientation on the elastic response for element 295 for 
a secondary orientation angle of 0=45® under 
mechanical loading only. The variations of the individual 
and normalized stress components with the primary 
orientation angle are shown in Figs. 5(a) and 5(b) 
respectively. None of the stress components shown in 
Fig. 5(a) varied significantly as i|r changed from 0® to 
10®. The normalized components, Fig. 5(b), exhibited 
some variation with 4r, however, these variations are not 
substantial, because the magnitudes of the respective 

1 The graphics software used to draw these figures superimposes the 
data symbols with the last symbol on top. The rest of the data 
symbols, as shown in the key for each figure, are underneath this 
symbol. 
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individual stress components are rather small. Thus, it 
should be realized that the variations in the normalized 
stress components with i|r, while indicative of the 
influence of the primary orientation on the elastic 
response of the SC PWA 1480, alone are not sufficient 
to assess the importance of the primary orientation 
effects. The variations in the actual magnitudes of the 



ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


individual components of stress also must be considered 
to determine the quantitative significance of the 
orientation effects. The stresses generated in elements 
245 and 445 under mechanical loading for 0=45° were 
similar to those presented in Fig. 5. In addition, for 
element 295, similar results were observed for 0=0°, 
under mechanical loading. 



0 2 4 6 8 10 

ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(a) INDIVIDUAL STRESS COMPONENTS. 



ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(b) NORMALIZED STRESS COMPONENTS 

FIG. 5: INFLUENCE OF PRIMARY ORIENTATION 
ANGLE ON STRESS COMPONENTS UNDER 
IMPOSED MECHANICAL LOADING FOR 
ELEMENT 295. 


(a) INDIVIDUAL STRESS COMPONENTS 



ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(b) NORMALIZED STRESS COMPONENTS 
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FIG. 6: VARIATION OF STRESS COMPONENTS 
WITH PRIMARY ORIENTATION ANGLE UNDER 
THERMAL LOADS FOR ELEMENT 445 (0=0°). 
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ANGLE OF PRIMARY ORIENTATION, y (Degrees) ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(a) INDIVIDUAL STRESS COMPONENTS (a) INDIVIDUAL STRESS COMPONENTS 



ANGLE OF PRIMARY ORIENTATION, y (Degrees) ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(b) NORMALIZED STRESS COMPONENTS (b) NORMALIZED STRESS COMPONENTS 

FIG. 7: VARIATION OF STRESS COMPONENTS FIG. 8: VARIATION OF STRESS COMPONENTS 

WITH PRIMARY ORIENTATION ANGLE UNDER WITH PRIMARY ORIENTATION ANGLE UNDER 

THERMAL LOADS FOR ELEMENT 445 (6=45°). THERMAL LOADS FOR ELEMENT 295 (6=0®). 

Figures 6(a)-(b) and 7(a)-(b) present the variations of the primary orientation angle for secondary orientations 

the individual and normalized stress components with of 0° and 45® under thermal loads along the X- 
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direction for element 445. The stress components in 
Figs. 6(a) and 7(a) show small variations with ijr except 
for 0 W and <j_ which exhibited noticeable variations 
with +. The magnitudes of the o^, and o a components 
for 6=45° are 16 percent and 41 percent larger, 
respectively, than the corresponding magnitudes for 
6=0°. This clearly indicates that a secondary orientation 
angle of 45° is undesirable for any primary orientation 
angle, +, between 0° and 10°. The normalized stress 
components, Figures 6(b) and 7(b), show slight 
variations as f increases for both secondary orientations 
considered (0=0° and 0=45°). 

Figures 8(a) and 9(a) show the stress components 
under thermal loading along Z-axis for element 295. 
The Oyy component is lower for 0=45° than for 0=0°, 
and both change as changes from 0* to 10°. The 
normalized stress components, Figs. 8(b) and 9(b), 

( T xyW(Vo,e and ( V*y( Vo, e > have exhibited lar S e 
variations with However, these variations, which are 
due to orientation changes as well as loading conditions, 
are not significant because of the relatively small 
magnitudes of these stress components. 

The stress components for element 295 under 
combined thermal and mechanical loading conditions 
were observed to obey the principle of superposition of 
linear elastic loadings. For this case, the normalized 
stress components for 0=0® and 0=45° exhibited 
similar variation with i|r as those observed in the case of 
pure thermal loading along the Z-axis, Figs. 8 (b) and 
9 (b). For brevity, figures of the individual and 
normalized stress components for the combined loading 
conditions for element 295 are not included. 

Figure 10 shows the stress components in element 445 
for combined thermal and mechanical loading with 
thermal loading applied along the X-axis as the primary 
orientation changes from i|r =0° to + = 10° , for 0=0® and 
0=45®. The changes in the stress components with i|r are 
small, which leads to the conclusion that the primary 
orientation does not have a significant effect on the 
elastic response under these conditions. However, the 
magnitudes of and a a , for a given were 41 and 29 
percent larger, respectively, for 0=45® than the 
corresponding magnitudes for 0=0®. As mentioned 
earlier for the loading case involving X-directional 
thermal load, the secondary orientation of 0=45® 
generated large stresses and hence should be avoided 
for thermal loads involving X-directional gradients. 



0 2 4 6 8 10 

ANGLE OF PRIMARY ORIENTATION, 6 (Degrees) 

(a) INDIVIDUAL STRESS COMPONENTS 

© 

© 



(b) NORMALIZED STRESS COMPONENTS 

FIG. 9: VARIATION OF STRESS COMPONENTS 
WITH PRIMARY ORIENTATION ANGLE UNDER 
THERMAL LOADS FOR ELEMENT 295 (0=45®). 

The normalized components for the combined thermal 
(along the X-direction) and mechanical loading case 
varied in a similar manner with ijr as in the case of pure 
thermal loading along the X-direction, Figs. 6(b) and 
7(b). It was found that the thermal gradient through the 
thickness had produced much larger stresses than that 
along the length (primary orientation of [001]) of the 
plate. This is because the thermal gradient through the 
thickness of the plate is much larger than that along Z- 
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ANGLE OF PRIMARY ORIENTATION, y (Degrees) 


(a) INDIVIDUAL STRESS COMPONENTS 
AT 0=0° 



(b) INDIVIDUAL STRESS COMPONENTS 
AT 8=45* 

FIG. 10: VARIATION OF STRESS COMPONENTS 
WITH PRIMARY ORIENTATION ANGLE UNDER 
THERMAL AND MECHANICAL LOADS FOR 
ELEMENT 445. 


axis. In addition, Lee. et al. (1991) analyzed a single 
crystal blade and the temperature bounds that they used 
corresponded with the data adopted to run this study, 
specifically, thermal gradients and magnitudes of 


temperatures. 


CONCLUDING REMARKS 

A parametric study was conducted using elastic finite 
element stress analysis to determine the influence of 
primary orientation angle, at two different secondary 
orientation angles of 0° and 45 s , on the elastic behavior 
of [001]-oriented nickel base single crystal superalloy, 
PWA 1480 under mechanical, thermal, and combined 
thermal and mechanical loads. The primary orientation 
angle was varied from 0° to 10° in increments of 2° and 
the stresses developed within a single crystal PWA 1480 
square plate under different loading conditions were 
computed. The following conclusions were drawn from 
this study: 

(1) The type of loading (mechanical, thermal, or 
combined) imposed on the square plate dictated the 
variation in the individual stress components with the 
primary and the secondary orientation angles. Variation 
of primary and secondary orientation angles did not 
influence the stresses generated within the single crystal 
PWA 1480 plate under the purely mechanical loading 
condition. However, under thermal and combined 
thermal and mechanical loading conditions, influence of 
both primary and secondary orientations was observed. 

(2) Accurate evaluation of the primary and secondary 
orientation effects should include consideration of the 
magnitudes of the individual stress components in 
addition to the normalized stress components for each 
loading condition. 

(3) The influence of the primary orientation angle, 
when constrained between 0° and 10°, on the elastic 
stresses generated within the single crystal PWA 1480 is 
substantially lower than the influence of secondary 
orientation angle, which is not typically controlled. 
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